Full-scale simulation study of a generalized Zakharov model for the generation of 

topside ionospheric turbulence 
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We present a full-scale simulation study of a generalized Zakharov model for the generation of 
the topside electrostatic turbulence due to the parametric instability during ionospheric heating 
experiments near the F region peak. The nonlinear tunneling of electromagnetic waves through the 
ionospheric layer is attributed to multiple-stage parametric decay and mode-conversion processes. At 
the bottomside of the F region, electrostatic turbulence excited by the parametric instability results 
in the conversion of the ordinary (O mode) wave into a large amplitude extraordinary (Z mode) wave 
tunneling through the F peak. At the topside interaction region, the Z mode undergoes parametric 
decay cascade process that results in the generation of the topside electrostatic turbulence and then 
conversion process yielding O waves that escape the plasma. This study may explain the observed 
topside ionospheric turbulence during ground based ionospheric heating experiments. 



I. INTRODUCTION 



There have been several radar observations of en- 
hanced plasma and ion lines at the topside ionosphere 
during ground based heating experiments in which the 
transmitting frequency has been kept somewhat lower 
than the critical plasma frequency of the F2 layer max- 
imum. A proposed explanation for the observed top- 
side plasma and ion lines is a nonlinear decay of large 
amplitude Z mode waves, which have been mode con- 
verted from the O mode heater wave at the bottom- 
side of the ionosphere, into Langmuir and ion acous- 
tic (lA) waves. The observations have attracted much 
interest since both the ordinary (O) and fast extraor- 
dinary (X) electromagnetic waves are reflected at or 
below the critical region where V = 1, except at a 
very narrow angular range around the critical angle of 
incidence (the Spitze angle) where the O mod e wave 

waves BuddeA. 1961 



mode converted to Z mode . _ _ 

Ginzburd . Il967/1970t \Gondarenko et all 120031: iMMfew, 



1990l |. fThroughout the manuscript, we use the standard 



notation V = oj^ /w^, Y = u Jr.pJio and U = Y"^ [Budderi . 
119611 : \Ginzbur<L Il967/1970j where LUpe is the electron 
plasma frequency, LUce is the electron gyrofrequency and lu 
is the transmitted frequency.) Enhanced topside plasma 
lines in the F2 region were first observed in heating exper- 
iments at Arecibo hy JGanqulv and Gordon 1983] . Dur- 
ing experiments where the EISCAT UHF beam was di- 
rected near the Spitze angle, i/sfeam et al\ 1990j | observed 



enhancements of ion lines in the topside interaction F2 
region where V = 1. Simil ar observation s at E ISCAT 
were done in the E region hy \Rietveld et al\ [2002j |. where 
there were no clear dependence of the power of the top- 
side echoes on the power of the plasma turbulence at the 
bottomside. There are several observations at EISCAT 
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\lsham et all 


1999|, 



In the experiment at EISCAT by 



many degrees from the Spitze angle, including vertical 
incidence. In order for normally incident O wave to tun- 
nel through the F peak, conversion of the O mode to the 
Z mode should occur, and the following conditions are 
necessary. First, that the transmitter frequency should 
be near the maximum plasma frequency of the F peak 
so that V < 1 + Y for t he peak plasma frequ ency of the 
F layer, as discussed by \Mishin et all 1997| | and shown 
in the linear E-region simulat i ons fo r a parabolic den- 
sity profile {Condarenko et all . l2003j |. Second, the ex- 
istence of the density /electric field fluctuations within 
the critical layer, which can occur due to the gener- 
ation of electrostatic turbulence excited by the para- 
metric instability or in the earlier p hase of the nonlin- 
ear t he rmal self-focu s ing in stability [Gondarenko et ail . 
l2006l |. \Mishin et ail (200 11 proposed that the topside 
plasma and ion lines could be due to resonant scatter 
of the HF pump wave into Z mode waves on the small- 
scale magnetic field aligned striations, which in turn are 
created by the HF pump wave via the resonant insta- 
bility, or, in the case of the low duty cycle experiments 
with pump frequencies near a multiple of the electron 
gyro frequency, due to a sec ond-order four wave interac- 
tion [Hua ng and j^Mcl . fl994 |. Ro cket observations du ring 
heating experiments in Arecibo [Gelinas et all . l2003l | re- 
vealed Z mode waves and field aligned striations above 
the O mode reflection point, and evidence of genera- 
tion of Bernstein waves where the local plasma frequency 
matches an electron gyro harmonic. 



During experiments where the EISCAT UHF beam 
was d irected along the magnetic field lines, \lsham et all 
1996j | discovered spectra with an unusually large spec- 



cycle (200 ms every 10 s), topside excitation of ion and 
plasma lines were also observed in a angular range of 



tral width (up to 200 kHz) and a mean frequency shifted 
by 200-300 kHz away from the heating frequency in both 
the upshifted and downshifted channels. This feature 
has been named the HF- induced outshifted line (HFOL). 
The HFOL spectra come from a height about 3-6 km 
above the refiection height. .Mishin et all |l997l | propose 
that the HFOL feature can be consistently explained if 
the upper region is the refiection layer of the Z mode ap- 
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pearing due to the mode conversion at the reflection level 
of the injected O mode wave. Their scenario includes 
large-scale density depletions at the Airy maxima of the 
Z mode, the parametric decay of t he Z mode into ele ctron 
Bernstein and lower hybrid waves \Kuo et al\ . \l99§^ . and 
the production of suprathermal electrons by the electro- 
static turbulence, which in turn excite Langmuir waves 
inside the density depletions, yielding the HFOL. 

The interaction between the electromagnetic pump 
wave and the electrostatic turbulence ha ve been 
modeled with _ Zakharov-type equations \Kuznetsoi\ . 



wave equation 



I1974I: \Zakharoi\ . 
\Fejer and KuoL 



II972II and wave k i netic equations 
[Feier and KuoV 1X9731 : \Perkins et dl . Il974{ . to un- 
derstand the wave spectra and radiation generation 
from electromagnetically driven Langmuir turbulence 
' DuBois etldl. Il99ll: [Caleev et all Il977l: \Mj0lhus et'oH . 



I1995I : \Stubbe et al\ . Il992l | . In this Letter we analyze the 
conversion of O mode waves into Z mode waves and the 
generation of topside ionospheric turbulence numerically 
with a full-scale simulations of a generalized Zakharov 
model [Eliasson and Stenfld . \200^ based on a nested grid 
method [Elias son. 200^^ We are using realistic length 
scales and ionospheric parameters relevant for the high- 
latitude facility EISCAT near Troms0, Norway, but sim- 
ilar parameters may be also found at high-latitude facili- 
ties such as HAARP, Alaska and HIPAS, near Fairbanks, 
Alaska. 



II. MATHEMATICAL MODEL 

In our one-dimensional simulation geometry, we as- 
sume that a large amplitude electromagnetic wave is 
injected vertically, along the z axis, into the vertically 
stratified ionospheric layer. The mathematical model, 
derived in Eliasso n and Sten flo [2007||, is based on a sep- 
aration of time scales, where the time-dependent quan- 
tities (here denoted tp{z,t)) are separated into a slow 
and a high-frequency time scale, ip = ips + iph- The 
high-frequency component is assumed to take the form 
iph = (l/2)-(/'(z, t) exp(— iu;oi)-l-complex conjugate, where 

represents the slowly varying complex envelope of the 
high-frequency field, and ojq is the frequency of the trans- 
mitted electromagnetic wave. Hence the time derivatives 
on the fast timescale will be transformed as d/dt 
d/dt — iujQ when going from the iph variable to the ip 
variable. The dynamics of the high-frequency electro- 
magnetic waves and electron dynamics is coupled non- 
lincarly to the slow timescale electron dynamics via the 
ponderomotive force acting on the electrons, and is me- 
diated to the ions via the electrostatic field on the slow 
timescale. 

We first present the envelope equations on the fast time 
scale. The transverse (to the z axis) components of the 
electromagnetic field in the presence of a slowly varying 
electron density is governed by the electromagnetic 
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dz^ 



£0 



(1) 



(2) 



where Aj_ and E_l is the transverse vector potential (in 
Coulomb gauge) and electric field, respectively. Here c 
is the speed of light in vacuum and Eq is the electric 
vacuum permittivity. The z component of the electric 
field is governed by the Maxwell equation 



dt 



£0 



(3) 



The high-frequency electron dynamics is governed by the 
continuity and momentum equation 



and 
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E_L + Ve X B( 



(5) 



no 



dz 



respectively, where Bq is the geomagnetic field, Ve is 
the electron collision frequency, e is the magnitude of 
the electron charge, rrie is the electron mass, v^e = 
the electron thermal speed, Ub is Boltz- 
mann's constant, is the electron temperature, and no 
is the background electron number density. We have as- 
sumed that the ions are immobile on the fast time scale 
due to their large mass compared to that of the electrons. 

The slowly varying electron number density is sepa- 
rated as = neQ{z) + ns{z,t) where neoiz) = nio{z) 
represents the large-scale electron density profile due to 
ionization/recombination in the ionosphere, and is the 
slowly varying electron density fluctuations. We regard 
the electrons as inertialess on the slow time scale, where 
we also assume quasi-neutrality rijs — n^s = ng. The dy- 
namics of the ions is governed by the continuity equation 



~dt 



no- 



dz 



0, 



(6) 



and, since the ions are assumed to be unmagnetized, the 
slow timescale plasma velocity is obtained from the ion 
momentum equation driven by the ponderomotive force 



dv,z 
dt 



Cldn^ 
no dz 



So 9|Ep 
Amino dz 



(7) 



where = [ksiTe + ST,)/™^]!/^ is the lA speed, is 
the ion temperature, mi is the ion mass, and we have 



denoted |E|2 = jE^p + \E,\'' and jE^p = \E,\'' ■ 
In obtaining Eq. ([7]), we have used the ion momentum 
equation 



mi 



3wf,j dui 



ar mi no oz 

together with the electron force balance equation 







ksTedns So 



hq dz Auq dz 



(8) 



(9) 



where we have neglected effects of the geomagnetic field 
Bq on the ponderomotive force (the |Ep-term) since we 
will consider high-frequency waves with frequencies much 
higher than the electron gyro-frequency. 




Ej^l (V/m) 



III. NUMERICAL SETUP 

We next define the simulation setup and physical pa- 
rameters used in the simulation; the details of the nu- 
merical implementat i on is described in the Appendix of 
\Eliasson and StenHd [2007 1. The simulation code uses a 
one-dimensional geometry, along the z axis. Our simu- 
lation box starts jit an altitude of 200 km and ends at 
400 km. While A_l and E_l are represented on a grid 
with grid size 2 m everywhere, the rest of the quanti- 
ties ne, Ve, rig. Viz and Ez are resolved with a much 
denser grid of grid size 4 cm at the bottomside and 
topside interaction regions regions, z = 286.2-287.2 km 
and z = 312.8-313.8 km, in order to resolve small- 
scale structures due to electrostatic turbulence. Out- 
side these regions, all quantities are resolved on the 
2-meter grid. This nested grid procedure is used to 
avoid a severe Courant-Friedrich-Levy (CFL) condition 
At < Sz/c on the timestep if the electromagnetic field 
is resolved on the dense grid 5z. The ionospheric den- 
sity profile will be assumed to have a Gaussian shape of 
the form nio{z) = no^maa; exp[-(z - Zmaa:)^/10^], where 

ill 



no.max — 5 X 10 m and z^ax = 300 x 10 m are 



the maximum density and the altitude of the F peak, 
and L = 31.6 x lO'^m is the density scale length. The 
external magnetic field Bq is set to 4.8 x 10"^ T and 
is tilted Q = 0.2269 rad (13°) to the vertical (z) axis 
so that Bo — Bo[xsm{9) — zcos(6')], which is the case 
at EISCAT in Troms0. We assume that the transmit- 
ter on ground has been switched on at t = 0. The 
simulation starts aX t = 0.67 ms, when the electromag- 
netic wave has reached the altitude 200 km. Initially, 
all time-dependent fields are set to zero, and random 
density fluctuations of order 10^ m"'^ are added to 
to seed the parametric instability in the plasma. The 
electromagnetic wave is injected from the bottomside of 
the simulation box at 200 km, by setting the x com- 
ponent of the electri c field to 1 V/m oii the upward 
propagating field (see lEliasson and Stenfld 2007l |). The 



transmitter frequency is set to loq = 3.66 x 10 s 
which is somewhat lower than the maximum F peak 



FIG. 1: The amplitude of the perpendicular (to the z axis) 
electric field in log scale and the density profile of the F layer 
(solid line). Electrostatic turbulence generated due to the 
parametric instability at the bottom and top sides of the F 
layer: at t ~ 1.7 ms, the O wave is transferred into Z waves 
which tunnel through the layer, and at t ~ 3.5 ms, the inverse 
conversion process yields O waves escaping the plasma. 



plasma frequency ujpe^max 
with no = no,max- On the other hand, the chosen 
wave frequency is higher than the maximum Z mode cut- 



off frequency ujz,max 
Y with V = 



3.58 X 10 s calculated from 



^„/^lmax, 



Y 



,r, 



eBo/nie ~ 8.5 x 10° s . By this choice of 



V =1 
and ujc 

pump frequency, the O mode is reflected by the over- 
dense plasma layer, while Z mode waves are allowed to 
propagate to the topside when electrostatic turbulence 
is generated due to the parametric instability. We use 
oxygen ions so that = 16 mp where irip = 1836 rrie 
is the proton mass. Further, = Ti = 1500 K so that 
VTe = 1-5 X 10^ m/s and Cg = 1.75 x 10^ m/s, and we use 
no ~ 3.3 X 10^^ m~^, which is the electron number den- 
sity a,t V — 1. The electron collision frequency is set to 
i^e = 10'^ s~^ and the effective ion "collision frequency," 
due to Landau damping, is set to = 2 x 10'^ s~^, where 
we have neglected that the ion Landau damping depends 
on the wavelength of the lA wave. The value used here 
is approximately valid for lA waves with wavelengths of 
-1-2 m. 



IV. NUMERICAL RESULTS 

The time evolution of the electric fields and ion den- 
sity fluctuations are shown in Figs. 1 and 2 from t = Q 
to t = 4.67 ms. Figure 1 shows the spatial and temporal 
evolution of the perpendicular (to the z axis) electric field 

After approximately 



\Ey\'- 



amphtude |E_l| = y\Ex 
1 ms, the electromagnetic wave reaches the bottomside 
interaction for the O mode at z « 287 km and for the 
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FIG. 2: Electrostatic turbulence generated due to the parametric instability shown with the electric field amplitude at (b) 
the bottomside and (a) topside interaction regions of the F layer at t ~ 1.7 ms and t ~ 3.5 ms, respectively. The ion density 
fluctuation at the interaction regions: (d) the bottomside when large-scale fluctuations are excited due to wave collapse and 
caviton formation at t ~ 1.7 ms and (c) the topside at t ~ 3.5 ms. 



X mode at a slighter lower altitude z « 270 km, after 
which a steady-state standing pattern quickly builds up. 
At t « 1.6 ms, there is an efficient conversion of the inci- 
dent O wave into the Z mode due to rapid development 
of the density fluctuations (discussed in connection to 
Fig. 2 below). The Z mode waves propagate from the 
bottomside sX z — 270 km to the topside interaction re- 
gion for the O mode at z — 312 km. At i = 3.5 ms, we 
observe excitation of O mode radiation at the topside at 
z w 312 km, and these O mode waves propagate away 
from the ionospheric profile on the topside. 

In Fig. 2, we have visualized the component of the 
high frequency electric field and the ion density fluctua- 
tions ns , in the vicinity of the topside and bottomside in- 
teraction regions for the O mode. In Fig. 2(b), we see the 
arrival of the O mode at the bott om side at ^ pa 1 ms, af ter 
which a standing wave pattern [Ginzburj}, '1967/1970] is 
set up. At t ~ 1.7 ms, the standing wave pattern breaks 
up into small-scale, large amplitude waves. At the same 
time, seen in Fig. 2(d), large amplitude ion density fluc- 
tuations are excited first in the region of the maxima of 
at z = 286.75km and z = 286.85km (see Fig. 2), and 
later spreading to lower altitudes. This transition into 
wave turbulence is closely correlated with the excitation 
of Z mode waves observed in Fig. 1. At i = 2-2.5 ms, 
we see in Fig. 2(a) the arrival of the Z mode at the top- 
side ionosphere. While some weaker waves propagate to 
higher altitudes, a band of large- amplitude electric field is 
set up at z « 313.2 k m where the Z mode wave meets the 
transformation point Mishin et all . l200ll : \Mj0lhui . Il990j | 
and converts into electrostatic waves. At t later than 
3.5 ms the Z mode wave decays into short wavelength 
Langmuir and lA waves due to a parametric instabil- 
ity, followed by Langmuir wave collapse and ion caviton 
formation. In Fig 2(c), the ion fluctuations are seen at 
the topside in a narrow band at 313.2 km. The large 
amplitude collapsed Langmuir waves and the associated 



ion cavitons have sizes of a few tens of centimeters or 
less and a spacing of a meter or le ss, similar as in the 
study of [M lasson ana Stenfld [2007l |. The topside Lang- 



muir turbulence gives rise to the radiation of O mode 
waves see n in Fig 1 at t > 3.5 ms, in a process described 
by \Mj0lhus et a/.l [l995( . 



V. DISCUSSION 

We have performed a simulation study of enhanced 
topside turbulence, where the transmitter frequency is 
kept slightly lower than the critical frequency of the F2 
layer. Our simulation study reveals that the electromag- 
netic wave decays parametrically on the bottomside of 
the ionosphere, in which collapsing wavepackets and ion 
density depletions are created. This is a three- wave decay 
of the O mode into one counter-propagating Langmuir 
(L) wave and one forward scattered lA wave, obeying the 
matching conditions uq = i-ol + ^lA of the frequencies 
and ko = ki, + kjA of the wavenumbers, where uJoiko), 
^^Likh) and uiAikiA) are given by the O mode, magne- 
tized Langmuir, lA dispersion relation, respectively. A 
wavelength of the O mode of the order 100 meters leads 
to wavelen gths of the Langmuir and lA waves of the order 
one meter \Eliasson and The linear insta- 



bility is followed by a rapid Langmuir waves collapse and 
the generation of ion cavitons. In the nonlinear phase 
of the instability, there is a strong generation of Z mode 
waves in the bottomside interaction region. These large- 
amplitude Z-mode waves are allowed to propagate into 
the denser plasma to the topside interaction region of the 
ionospheric layer, where they reach their transformation 
point and their amplitude increase as their group speed 
decreases and they are turning into electrostatic waves. 
At the topside, the Z mode waves undergo a paramet- 
ric instability into small-scale (few tens of centimeters) 
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coupled electrostatic and lA waves, followed by Lang- 
muir wave collapse and the formation of ion cavitons. In 
this la st step, Langinuir tu rbulence generated O mode 
waves [Mj0lh us et ai\ . Il995l | are propagating away from 
the plasma layer on the top side. We note that the wave 
tunneling and the excitation of topside turbulence occur 
within a few milliseconds, and may explain the observed 
topside enhanced plasma and ion lines far away from the 
Spitze region, in particular during low duty-cycle iono- 
spheric interaction experiments where the impact of iono- 



spheric irregularities are minimized \lsham et all . Il999l | . 
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